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Abstract 

An infinite dimensional àlgebra, which is useful for deriving exact solutions of the gen- 
eralized pairing problem, is introduced. A formalism for diagonalizing the corresponding 
Hamiltonian is also proposed. The theory is illustrated with some numerical examples. 
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Pairing has long been considered an important interaction in physics. The concept 
can be traced back to the seniority scheme introduced by Racah in atòmic physics. 1 Its 
physical significance was first realized in the study of superconductivity. 2 Following the 
suggestions of Bohr, Mottelson, and Pines, 3 the first detailed application of pairing in 
nuclei was made by Belyaev. 4 The concept has since been applied to other phenomena: 
high T c superconductivity, 5 ' 6 applications using the Hubbard model, 7 pairing phenomena 
in liquids, 8 and metal clusters. 9 

BCS methods have yielded major successes in studies of superconductivity. When 
applied to nuclei, however, some negatives come with the positives. First of all, not 
only is the number of nucleons in a nucleus typically small, the number of valence par- 
tides (n ~ 10) which dominates the behaviour of low-lying states is too few to support 
underlying assumptions of the BCS approximation, specifically, 5n/n is not negligible. 
As a consequence, particle-number-nonconservation effects enter and can lead to serious 
difliculties, such as spurious states, nonorthogonal solutions, etc. Secondly, an essential 
feature of pairing correlations are even-odd differences, which are driven mainly by Pauli 
blocking. It is difncult to treat these differences in the BCS formalism because different 
quasi-particle bases must be introduced for different blocked levels. For these reasons, a 
particle-number-conserving method for treating the pairing problem has been suggested 
for well-deformed nuclei. 10 The method uses a configuration energy truncation scheme, 
and takes the strength of the pairing interaction to be the same for all orbitals. In this 
limit the pairing Hamiltonian can be diagonalized in truncated configuration spaces. Be- 
cause the theory applies to well-deformed nuclei, each orbital can only be occupied by a 
single pair of partides. 

The generalized pairing Hamiltonian for spherical nuclei can be written as 

H = J2e 3 a] m a rrn - \G\ £ c jr S + (j)S-(j), (1) 

jm jj' 

where the €j are single-particle energies and S ± (j) and S°(j) are the pairing operators for 
a single-j shell defined by 

m>0 

S~(j) = £ (-r m a,_ m a jm , 

m>0 

S °Ü) = \ Yl ( a ]m a jm + a)_ m a ó _ m - 1) = )-(Nj - fy-), (2) 

Z m>0 Z 

where Qj = j + 1/2 is the maximum number of pairs in the j-th shell, Nj is the particle 
number operator for the j-th shell, and the Cjj/ measure the pairing strength between 
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different j-shells. In general, for N pairs, Hamiltonian (1) can be diagonalized in bases 
states that are products of the single-j shell pairing wave functions: 

\N >=Y.B klk2 ... kp S+ k >S+ k * ..-SÍ fc l0 >, (3) 

where the summation is restricted by 

í> = ^ (4) 
i=\ 

the B klk2 ... kp are expansion coefficient that need to be determined, and |0 > is the pairing 
vacuum state which satisfies the condition 



Sj\Q> =0 for all j. (5) 



The dimensionality of the Hamiltonian matrix in this basis increases very rapidly with 
increasing iV and the number of shells p. It is less than or equal to the dimension of 
the irreducible representation (irrep) [NÒ] of the unitary group U(p) due to the Pauli 
principle, 

dim < tE+tLzM. ( 6 ) 

~ N\(p-1)\ V ; 

The equal sign holds in (6) when all the single-j shell pairing wave functions in the 
summation of (3) are Pauli allowed. Indeed, the problem quickly becomes intractable 
because there are no analytical expressions or recursion relations for determining the 
B kl k 2 -k p coeficients. 

Following the quasi-spin approximation, 11 consider a simpler Hamiltonian 

H = -\G\S+S», (7) 

where 

So + = E ^ + U), S = £ Cj S-{j) (8) 
j j 

with the coefficient q, c* satisfying the condition 



3 



EM 2 = i- (9) 

i 

This defines generalized pairing as proposed by Talmi. 12 Clearly, in the notation of (1), 
Cjji = c*Cjf with \cj\ 2 giving the percentage of single-j shell pairing in the Hamiltonian. 
In what follows, the Cj are taken to be real. 

To diagonalize Hamiltonian (7), consider an àlgebra generated by 

^ = Es 2m+1 ^0')- ( 10 ) 
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It is easy to show that these generators satisfy the following commutation relations: 



m+n+l i 

[S® n , Si^] = ± 'S'm+n· (H) 

Therefore, the {S^, [i = 0, +, — ; m — 0, ±1, ±2, • • ■} form an infinite-dimensional àlgebra, 
one that differs only slightly from a general Lie àlgebra of the affine type. 

The unique lowest-weight state of this àlgebra is simply the product of the single-j 
shell pairing vacua with arbitrary seniority quantum numbers. Therefore, it suffices to 
consider the total seniority zero case. The lowest-weight state satisfies 

,S m |0>=0; m = 0, ±1, ±2, (12) 

and 

s° m \o >= —Ei^r^-i >= A -i° > • ( i3 )- 

Z 3 

Furthermore, it can be proven that the eigenvectors of H for any N and non-zero energy 
eigenvalue can be expanded as 

\N >= ^ xTxT ■ ■ ■ xTtS+S+tS* ■ ■ ■ S: N J0 >, (14) 
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where 



rii 



-1, - 2, • • • if c)xí > 1, 
0, 1, 2, ••• if cja* < 1. 



In any case, up to a normalization constant (14) can always be written as 



where 



\N >— SqS^S^ 2 ■ ■ ■ S'^ v _ 1 |0 >, 



(15) 



(16) 
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While xi, x 2 , ■ ■ ■ ,%n-i are real or complex numbers satisfying the relation 



(17) 



^Sfycfo - 1 2 =~ + E^^> i = l, 2,.-., iV-l, 



1 " 



with 



where 



JV-l j 



yi = Xi/a, a 



h + 2A 1 



, h = S/(-|G|). 



(18) 



(19) 



(20) 



Therefore, the coemcients Xj (i=l, 2,- • •, N-l) and eigenvalues of the pairing energy E ^ 
are simultaneously determined by the system of equations (18) and (19). 



While eigenvectors for E = can be expanded as follows 

1^ >= E X T X 2 2 ' ' ' X 'n' ^n^m ' ' ' ^n N \^ >> 



(21) 



Ui 



where the restrictions on the integers rii are the same as those given by (15), and the 
expansion coemcients Xi are determined by the following set of equations 



E^— 



x x c) 



for N = 1, 



(22) 
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c? _ i 



2íí jT _i-^ = 2^— , i = 1, 2, • • • , iV, for N >2. (23) 
It is instructive to write down the first few energy eigenvalues and eigenstates. 

fc|0 >= 0; (24) 

h\i>=j2\c j \ 2 tt j \i>, |i >= s+io >, 

3 

h\N = l, p> = 0, \N = l, p>=^ _ C \ 2 S+\0>, p = 1, 2, (25) 

where x p is determined by (22); 

fe|2>= /i|2>, |2 >= S'+S'+|0 >, (26) 
where h is one of the solutions of the equation 

l t^ 3 x-^— = l (27) 



where 



2£í J J xc|-l 



*=-^v (28) 

h|7V = 2, p > = 0, p = 1, 2, ••• (29) 



where 



|7V = 2, p >= J2 T , v^A VY\ S t S t^ >» ( 3 °) 

fjí (l-xícj)(l-x^cj,) J J 

and xf are determined by (23). 



Energy levels of the generalized pairing interaction for the j = 1/2, 3/2, and 5/2 
case with ci/2 = v^TÏ, c 3 / 2 = a/cT3, and c 5 / 2 = x/ÜlJ are shown in Fig. 1. It can be seen 
from Table 1 that the energy level with the largest Cj value for the highest j orbit is the 
lowest one for any fixed N. Finally, an example of the lowest levels for given N in the 
5-th and 6-th shell, respectively, are shown in Table 2. 



6 



Acknowledgment 



Supported by the National Science Foundation through Cooperative Agreement No. 
EPS-9550481 and Grant No. 96030006. 



[I] G. Racah, Phys. Rev. 62, 438 (1942) 

[2] J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. Rev. 108, 1175 (1957) 

[3] A. Bohr, B. R. Mottelson, and D. Pines, Phys. Rev. 110, 936 (1958) 

[4] S. T. Belyaev, Mat. Fys. Medd. 31, 11 (1959) 

[5] M. Randeria, J. M. Duan, and L. Y. Shieh, Phys. Rev. Lett. 62, 981 (1989) 

[6] S. Schmitt-Rink, C. M. Verma, and A. E. Ruckenstein, Phys. Rev. Lett. 63, 445 (1989) 

[7] C. N. Yang, Phys. Rev. Lett. 63, 2144 (1989) 

[8] D. W. Cooper, J. S. Batchelder, and M. A. Taubenblatt, J. Coll. Int. Sci. 144, 201 (1991) 

[9] M. Barranco, S. Hernandez, and R. J. Lombard, Z. Phys. D22, 659(1992) 

[10] J. Y. Zeng and T. S. Cheng, Nucl. Phys. A405, 1 (1983); A411, 49 (1984); A414, 253 (1984) 

[II] A. K. Kerman, Ann. Phys. (N. Y.) 12, 300 (1961) 
[12] I. Talmi, Nucl. Phys. A172, 1 (1971) 



7 



Table 1 . The lowest energy levels for j = 

Case a: c 1/2 = \/0T, c 3/2 = 
Case b: c 1/2 = V0~2, c 3 / 2 = 
Case c: cy 2 = c 3 / 2 = 

Case d: c 1/2 = Vol, c 3/2 = 



1/2, 3/2, 5/2 with: 

c 5/2 = Vol; 
/OT, c 5/2 = v/oTr; 
/ 0l, c 5/2 = v/oT; 

/02, c 5/2 = 70l. 



/» = S/(-|G|) 



N 


a 


b 


c 


d 

















1 


2.6 


2.5 


1.9 


1.4 


2 


4.083 


3.866 


2.858 


2.136 


3 


4.735 


4.441 


3.278 


2.481 


4 


4.674 


4.377 


3.100 


2.320 


5 


4.084 


3.866 


2.858 


2.136 


6 


2.600 


2.500 


1.900 


1.400 



Table 2. Lowest energy levels for 5th and 6th shells. The parameters are chosen as follows: 
Ci/2 = V0.2, c 3 /2 = \/0.3, c 5 / 2 = \/0.22, Cg/ 2 = \/0.28 for 5th shell, whereas Ci/ 2 = V0.2, 
c 3 / 2 = \/0.14, c 5 / 2 = \/0.22, c 7 /2 = y/0.28, c n / 2 = \/0.l6 for 6th shell, respectively. 



(N, h = E/{-\G\)) 



5th shell (0, 0) 

(6, 9.292) 



6th shell 



(0,0) 

(6, 13.035) 
(12, 11.655) 



(1, 2.860) 
(7, 9.024) 

(1, 3.220) 
(7, 13.778) 
(13, 10.096) 



(2, 5.192) 
(8, 8.243) 

(2, 6.014) 
(8, 14.125) 
(14, 8.167) 



(3, 6.999) 
(9, 6.954) 

(3, 8.387) 
(9, 14.082) 
(15, 5.874) 



(4, 8.283) 
(10, 5.515) 



(5, 9.046) 
(11, 2.86) 



(4, 10.345) (5, 11.893) 
(10, 13.653) (11, 12.843) 
(16, 3.220) 
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N = 1 


N = 2 


N = 3 



Fig. 1. An example of excited energy levels for the generalized pairing interaction with N < 3, 
where the number on the left indicates the degeneracy when it is greater than one. 



